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摘要 : 使 用 不 同 于 传统 的 类 全 里 叶 方法 的 celerite FIAT HTH EK AC 01.02 的 y 射线 和 射 
电波 段 的 长 期 光 变 特征 。 结 果 表 明 celerite 中 最 简单 的 阻尼 随机 游 走 模型 (damped random 
walk, DRW) 可 以 成 功 地 拟 合 射电 和 7 射线 的 长 期 光 变 曲线 ， 而 较为 复杂 的 二 阶 随 机 过 程 
(stochastically-driven, damped harmonic oscillator, SHO) 并 没有 显著 提升 拟 合 优 度 。 得 到 的 内 
ZEB 7 射线 光 变 特征 时 标 约 为 3 a， 如 此 长 的 时 标 不 能 在 轻 子 辐射 模型 中 产生 ， 但 是 在 强 子 模型 
中 是 可 能 存在 的 ， 由 此 推测 4C 01.02 长 期 的 y 射线 辐射 可 能 起 源 于 强 子 过 程 。 射 电光 变 的 内 豪 
特征 时 标 约 为 10 a， 它 可 能 对 应 于 射电 辐射 区 的 逃逸 时 标 ， 这 表明 射电 辐射 产生 于 大 尺度 喷 流 。 
X 键 gu: 活动 星系 核 ， 兆 变 体 ， 射 电光 变 ， 射线 光 变 ，celerite 方法 
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(BL Lacertae object, BL Lacs) 和 平 谱 射 电 类 星体 (flat spectrum radio quasars, FSRQs). 388 
常 ， 蝎 虎 座 型 天 体 本 质 上 是 低 喷 流 功 率 的 射电 星系 ， 而 平 谱 射电 类 星体 本 质 上 是 高 喷 流 功率 
的 射电 星系 ””。 洽 变 体 的 辐射 履 盖 了 整个 电磁 频谱 ， 主 要 由 相对 论 性 喷 流 的 辐射 主导 。 炊 
变 体 典 型 的 多 波段 能 谱 分 布 有 两 个 不 同 的 峰 ， 低 能 量 峰 在 红外 到 XX 射线 波段 ， 高 能 量 峰 在 
7 射线 能 量 段 ”。 沪 变 体 是 最 重要 的 河 外 y 射线 源 。 
WEE PR 7 射线 辐射 可 由 轻 子 过 程 和 强 子 过 程 产 生 ， 目 前 没有 定论 。 在 轻 子 喷 流 模型 
H, y 射线 由 喷 流 中 相对 论 性 电子 对 低能 光子 的 逆 康 普 顿 散射 产生 ， 包 括 同步 自 康 普 顿 散 
射 ”“ 和 外 康 普 顿 散 射 "”。 在 强 子 模型 中 ，7 射线 主要 来 源 于 喷 流 中 高 能 质子 的 同步 辐射 或 
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者 质子 与 光子 相互 作用 产生 的 辐射 。 
次 变 体 的 辐射 有 明显 的 光 变 。 在 7 射线 波段 ， 光 变 时 标 可 以 短 至 几 分 钟 量 级 ”， 而 
在 射电 和 光学 波段 光 变 时 标 可 以 长 至 数 十 年 ”。 光 变 的 分 析 是 探索 活动 星系 核 物理 本 
质 的 重要 手段 ””。Fermi y 射线 空间 望远镜 已 经 运行 了 十 多 年 ， 搭 载 在 其 上 的 Large 
Area Telescope (LAT) 提供 了 明亮 炊 变 体 的 高 质量 的 长 期 光 变 曲线 。Owens Valley Radio 
Observatory (OVRO) 的 40 m 望远镜 对 光 变 体 的 射电 活动 进行 监测 ， 补 充 了 Fermi 的 7 射 
线 观 测 ”。 这 项 工作 的 目的 是 利用 OVRO 和 LAT 数据 研究 高 红 移 耀 变 体 AC 01.02 的 射电 
和 7? 射线 长 期 光 变 特性 。 
随机 过 程 模型 (也 称 为 高 斯 过 程 模型 ) 是 分 析 活 动 星 系 核 光 变 的 有 力 工具 ”。 活 动 
星系 核 吸 积 盘 的 光学 光 变 可 以 用 一 个 简单 的 随机 过 程 模型 描述 ” ， 即 阻尼 随机 游 走 模型 
DRW)。 最 近 ， 随 机 过 程 方法 被 用 于 分 析 泡 变 体 的 y 射线 光 变 ”™， 这 些 作 者 采用 了 由 
Kelly 等 人 ”开发 的 CARMA 模型 (continuous time autoregressive moving average), fti] 
声称 CARMA(1,0) 过 程 (BI DRW 模型 ) 不 能 很 好 地 描述 某 些 焰 变 体 的 7 射线 数据 的 特征 。 
= celerite 是 一 个 新 的 随机 过 程 模型 ”， 近 几 年 被 广泛 地 应 用 到 了 天 体 的 光 变 研究 中 。Yang 
2 a A 用 celerite 模型 分 析 了 27 个 可 能 存在 y 射线 准 周期 震荡 的 耀 变 体 的 LAT 数据 ， 仅 在 
其 中 两 个 耀 变 体 中 发 现 了 准 周期 震荡 可 能 存在 的 证 据 。Zhang “A 使 用 celerite 在 PKS 
0521-36 中 发 现 了 y 射线 准 周期 震荡 存在 的 显著 证 据 。 这 些 工 作 表 明 celerite 模型 适合 用 了 
耀 变 体 的 光 变 研究 。 
这 里 我 们 使 用 celerite 模型 来 模拟 炮 变 体 AC 01.02 的 射电 以 及 yy 射线 的 光 变 曲线 。 除 
了 DRW 模型 ， 我 们 还 采用 一 个 更 高 阶 的 阻尼 简 谐振 荡 模 型 (SHO)， 来 研究 SHO 模型 是 否 
能 比 DRW 模型 更 好 地 描述 光 变 。 此 外 我 们 也 探讨 了 射电 和 射线 光 变 的 物理 起 源 。 本 文 
第 2 章 简单 介绍 我 们 使 用 的 随机 过 程 模型 ， 第 3 章 详细 展示 分 析 得 到 的 拟 合 结果 ; 最 后 在 
第 4 章 对 得 到 的 结果 进行 讨论 并 给 出 总 结 。 


2 随机 过 程 模型 


= 随机 过 程 (也 称 为 高 斯 过 程 ) 是 由 一 个 均值 函数 pola) 和 一 个 协 方差 ， 或 者 叫 核 函数 
© Ka (Wry) 组 成 的 随机 模型 中 , 这 里 的 0 和 Q 是 参数 。 此 处 ， 坐标 X= (ea ae susce gt 
处 的 数据 y = (yy? y) 的 对 数 似 然 函 数 是 : 


1 1 N 
InL(0,a) = In p(y|X,0,0) = -370 Ka ro — sn detK, — >In (27) , (1) 


AF, ply|X,0,a) 是 似 然 函数 ，rg RANE Ka Æ N x ON 的 协 方差 矩阵 ， 其 中 N 为 
数据 个 数 。 

对 于 给 定 的 数据 (y, X); BBO 和 a 可 以 使 用 非 线 性 优化 方法 使 似 然 函 数 L(0, a) 最 
大 时 得 到 最 优 值 ，9 和 a 的 误差 可 以 通过 设置 一 个 先 验 的 p(0, a)， 而 后 通过 马尔 科 夫 链 蒙 特 
卡 洛 (Markov chain Monte Carlo, MCMO) 算法 计算 的 参数 后 验 分 布 得 到 。 
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DRW UTE: ERO s 
k(t) = ae-t/romw , (2) 
其 中 a 是 振幅 参 数 ， TDRW 是 此 模型 的 特征 时 标 。 
SHO 模型 的 微分 方程 为 : 
d? d 1 
[inq artem o 


其 中 ，wo 是 阻尼 振子 的 角 频 率 ， 可 以 转化 为 特征 时 标 只 Ho = 2r/w Q 是 振子 的 品质 因 
d. e(t) 则 是 白 噪声 ， 对 应 的 功率 谱 (power spectral density, PSD) 为 


2 Sows 
= emer sor | M 


RH, So 是 w = wo 时 的 功率 
这 里 ， 我 们 采用 celerite 模型 来 拟 合 光 变 | 


deviance information criterion) 进行 模型 选择 。 


线 ， 并 用 贝 叶 斯 信息 准则 (Bayesian criterion 


3 4C 01.02 光 变 曲线 拟 合 结 


4C 01.02 (PKS B01064-013) 是 一 个 高 红 移 (z = 2.099) 的 明亮 y 射线 平 谱 射电 类 上 
Wk. LAT 数据 的 分 析 遵循 点 源 分 析 的 标准 @@。 我 们 选用 大 于 100 MeV 的 数据 ， 用 标准 工具 
Fermipy 以 30 d 数据 为 一 个 点 来 产生 光 变 曲线 。TS < 292 的 数据 点 被 排除 ， 以 确保 在 分 
析 中 采用 的 数据 可 靠 。 

AC 01.02 的 15 GHz 射电 光 变 曲线 由 OVROB 提供 。> 射线 的 数据 覆盖 MJD 54697 到 
59377 (2008 年 8 月 19 HÆ 2021 ^E 6 A 12 H) 期间， 射电 数据 覆盖 MJD 54910 到 58071 
(2009 ££ 3 H 12 H $1 2017 4E 11 A 14 A) 期 间 。 

我 们 用 celerite 模型 中 的 DRW 和 SHO 模型 拟 合 光 变 曲线 ， 并 使 用 MCMC sampler 
emcee 来 提升 拟 合 质量 。 每 次 拟 合 我 们 让 MCMC sampler 使 用 32 个 平行 的 行走 器 游 走 
10000 步 作 为 预 拟 合 ， 再 用 20000 步 来 得 到 MCMGC 的 样本 。 
图 四 展示 了 y 射线 光 变 曲线 的 DRW 模拟 结果 ， 图 2 为 SHO 的 模拟 结果 。 对 于 每 个 
模型 ， 我 们 展示 了 光 变 曲线 的 拟 合 ， 标 准 化 残 差 的 自 相 关 函 数 (auto correlation function, 
ACF)， 标 准 化 残 差 分 布 以 及 由 拟 合 结果 计算 得 到 的 功率 谱 。 标 准 化 残 差分 布 中 的 红色 虚线 
是 对 标准 化 残 差 的 常数 最 优 拟 合 。 由 MCMC 得 到 的 模型 参数 的 后 验 分 布展 示 在 图 e)。 


BRI 


?http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Pass8 usage.html 


示 区 域 的 有 点 源 时 模型 的 最 大 似 然 函 数 ，Lo 则 是 没有 源 时 


?'Test Statistics (TS) 定义 为 2lg(L/Lo)， 其 中 工 是 对 在 目 书 
的 最 大 似 然 函数 中 
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ik: a) 光 变 曲线 最 优 拟 合 与 LAT 数据 ， 阴 影 区 域 为 95% 置信 区 间 ; b) 标准 化 残 差 自 相 关 函 数 (ACE): c) 标 
准 化 残 差 分 布 (红线 为 最 优 常数 拟 合 ); d) 功率 谱 ，e) corner 模型 的 参数 分 布 。 


图 1  celerite 4C 01.02 的 射线 光 变 曲线 的 DRW 拟 合 结果 
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Hs 各 图 含义 与 图 四 相同。 


2 celerite 4C 01.02 AY y 射线 光 变 曲线 的 SHO 拟 合 结 


从 图 中 可 以 看 出 ，DRW 和 SHO 的 拟 合 结果 都 是 可 以 接受 的 : 残 差 的 自 相关 函数 几乎 
都 在 白 噪声 95% 置信 区 间 的 阴影 范围 内 ， 说 明 两 种 拟 合 都 抓 住 了 光 变 曲线 中 的 特征 。 此 外 ， 
从 残 差分 布 图 中 可 以 看 到 残 差 几乎 均匀 随机 地 分 布 在 y = 0 的 两 人 出， 而 且 常 数 拟 合 的 结果 非 
常 接近 0 ( 约 化 卡 方 值 x? 2 1)， 这 也 表明 这 两 个 模型 的 拟 合 结果 都 是 可 以 接受 的 ，DRW 和 
SHO 的 模型 参数 都 限制 得 很 好 。 

图 3 和 图 4 展示 了 射电 数据 的 拟 合 结果 ， 与 7 射线 的 结果 相似 ， 两 个 模型 都 成 功 地 模 
拟 了 光 变 曲线 ， 模 型 参数 也 限制 得 很 好 。 
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iS 各 图 含义 与 图 四 相同 。 


3 celerite 4C 01.02 的 OVRO 射电 数据 的 DRW WAAR 
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注 ， 各 图 含义 与 图 四 相同 。 


图 4 celerite 4C 01.02 的 OVRO 射电 数据 的 SHO WAAR 


为 了 定量 区 别 DRW 与 SHO， 我 们 计算 两 个 模型 的 DIC. DIC 的 定义 如 下 : 


DIC = pp 4 D(0) . (5) 


tb, DO 是 方差 DO) 的 平均 值 ， 等 于 —218(L(0) pp 是 参数 的 有 效 数 ， 用 于 描述 模型 
的 复杂 程度 ， 我 们 采用 pp = 2 var(D(0).. DIC 随 拟 合 优 度 的 增加 而 减 小 ， 随 模型 复杂 度 
的 增加 而 增 大 。 
从 表 回 中 我 们 得 知 DRW 与 SHO 的 DIC 相差 很 小 ， 这 表明 这 两 种 模型 的 拟 合 结果 是 
可 比较 的 。 
在 表 加 中， 我们 列 出 了 两 种 模型 计算 得 到 的 特征 时 标 。 可 以 看 出 ， 对 于 y 射线 数据 ， 
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表 1 DRW 与 SHO 模型 的 DIC 
源 名 称  DICsno, DICprw,,  DlCsuo,  DlICpmw, 
4C 01.02 26 24 827 831 


TDRW,y 和 Tsuo,, TEVAZETREEL VI eB, y 射线 的 特征 时 标 大 概 在 200 到 300 ds MATA 
电波 段 ， TSHO,r 大 于 TDRW,re 射 电 数 据 覆 盖 的 时 间 略 长 于 TSHO? 这 可 能 会 在 确定 TSHO,r 时 
引入 较 大 的 偏差 ， 射 电 数据 的 长 度 大 概 是 rprw 的 3 倍 ， 所 以 rprw, 基本 上 是 可 信 的 ”。 
Tpnw., 大 概 在 1000 d EA, HERT y 射线 时 标 。 


表 2 最 优 拟 合 特征 时 标 d 
源 名 称 ”Ts8HO,y TDpRw TSHO,r TDRW,r 


4C 01.02 369135 2451420 24161363 1075+435 


4 讨论 与 总 结 


在 这 项 工作 中 ， 我 们 使 用 celerite PLAY IW A WEAR HK 4C 01.02 在 射电 和 y 射线 波段 的 
流量 变化 ， 并 得 到 每 个 波段 长 期 变化 的 特征 时 标 。 结 果 表 明 ， 最 简单 的 随机 模型 DRW 和 
较为 复杂 的 SHO 模型 都 能 成 功 描述 射电 和 射线 的 光 变 曲线 ， 两 者 的 拟 合 结果 相当 ， 因 此 
更 复杂 的 SHO 模型 没有 必要 使 用 。 这 与 活动 星系 核 吸 积 盘 光 变 的 行为 类 似 呈 。 值 得 注意 的 
是 ， 我 们 从 光 变 曲线 中 得 到 的 特征 时 标 是 观测 值 ， 内 豪 时 标 应 为 7? = r 6p/(L 十 2)， 其 中 
óp 是 多 普 勒 因子 。4C 01.02 的 多 普 勒 因子 是 未 知 的 ， 平 谱 射 电 类 星体 的 平均 多 普 勒 因子 为 
7~15”， 因 此 我 们 采用 平 谱 射 电 类 星体 的 典型 多 普 勒 因子 dp = 10， 修 正 的 y 射线 特征 
时 标 大 约 为 3 a， 而 射电 波段 的 特征 时 标 约 为 10 ac 
通常 认为 这 个 光 变 特征 时 标 对 应 于 辐射 系统 的 菜 个 时 标 ， 如 冷却 时 标 、 加 速 时 标 或 逃逸 
时 标 等 ”。 喷 流 中 的 y 射线 辐射 可 由 轻 子 过 程 和 强 子 过 程 产生 ， 目 前 没有 定论 “。4C 01.02 
是 一 个 典型 的 平 谱 射 电 类 星体 ， 在 轻 子 模型 中 ， 它 的 y 射线 通常 由 外 康 普 顿 过 程 产生 ， 对 
应 的 冷却 时 标 约 为 10 d 量 级 ”。 这 个 冷却 时 标 远 小 于 我 们 得 到 的 光 变 特征 时 标 ， 因 此 我 们 
认为 y 射线 光 变 特征 时 标 不 能 由 轻 子 辐射 效应 产生 。 而 在 强 子 模型 中 ， 相 关 过 程 的 冷却 时 
标 可 达 数 年 甚至 上 百年 ”。 这 表明 我 们 发 现 的 长 达 约 3 a 的 y 射线 特征 时 标 可 以 由 强 子 过 
程 产生 ， 即 4C 01.02 的 y 射线 辐射 可 能 起 源 于 强 子 过 程 。 具 体 的 光 变 起 源 需要 在 含 时 强 子 
模型 框架 下 继续 研究 。 

喷 流 中 的 射电 辐射 由 高 能 电子 的 同步 辐射 贡献 ， 同 步 辐射 过 程 的 冷却 时 标 主 要 由 磁场 
决定 。 在 典型 的 喷 流 磁场 的 参数 空间 中 ， 同 步 辐射 的 冷却 时 标 可 达 10 as ， 因 此 射电 的 特 
征 时 标 表 明 射 电光 变 可 能 起 源 于 辐射 区 磁场 的 扰动 。 此 外 ， 射 电 的 特征 时 标 也 可 能 对 应 于 辐 
射 区 的 逃逸 时 标 ， 由 此 可 以 限制 射电 辐射 区 的 半径 约 为 10 cm， 即 射电 辐射 产生 于 大 尺度 
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喷 流 ， 这 也 符合 一 般 的 认识 。 值 得 强调 的 是 对 oy 射线 辐射 起 源 的 推断 是 针对 长 期 光 变 特征 
的 ， 不 一 定 适用 于 短期 的 快速 光 变 。 
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Studying Radio and y-ray Variability Characteristic 
Timescales of the Blazar AC 01.02 with Stochastic Process 
Method 
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Abstract: Here, we use the celerite method, which is different from the traditional Fourier- 
like methods, to analyze the y-ray and radio long-term variabilities of 4C 01.02. The results 
show that the simplest kernel in celerite, damped random walk (DRW), can successfully fit 
both radio and y-ray long-term light curves, and the more complex second-order stochas- 
tic process (stochastically-driven damped harmonic oscillator, SHO) does not significantly 
improve the goodness of fit. 

The intrinsic characteristic timescale of y-ray variability is about 3 a, and such a long 
timescale cannot be generated in the leptonic emission model, but it is allowed in the hadronic 
model. We therefore speculate that the long-term y-ray emission of 4C 01.02 may originate 
from hadronic process. The intrinsic characteristic timescale of radio variability is about 10 
a. It may correspond to the escape timescale of the radio emission region, indicating that 


radio emission is produced in the large-scale jet. 


Key words: active galactic nucleus; blazar; radio variability; y-ray variability; celerite 
method 


